Noisy Kuramoto-Sivashinsky equation for an erosion model 
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We derive the continuum equation for a discrete model for ion sputtering. We follow an approach 
based on the master equation, and discuss how it can be truncated to a Fokker-Planck equation 
and mapped to a discrete Langevin equation. By taking the continuum limit, we arrive at the 
Kuramoto-Sivashinsky equation with a stochastic noise term. 
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I. INTRODUCTION 

The field of nonequilibrium interface growth is cur- 
rently being investigated by a variety of approaches . 
One direction of research is the formulation of Langevin 
equations which are studied by means of analytical and 
numerical methods. Typically, the terms appearing in 
Langevin equations are derived by the use of symmetry 
arguments. A different approach consists in the formula- 
tion of discrete models which contain the relevant phys- 
ical mechanisms present in the problem at hand. The 
models are studied by computer simulations in order to 
determine which universality class they belong to. One 
can try to relate the models to continuum descriptions. 
In many occasions such an identification is indirect since 
it is difficult to derive the continuum description for a 
given set of microscopic growth rules. 

There exist cases in which it is possible to derive 
Langevin equations describing various discrete growth 
models The approach followed in these cases is 

based on the master equation which determines the evo- 
lution of the joint probability density P{H,t), where H 
specifies the interface configuration and t denotes time. 
The master equation reads 



dP{H,t) 
di 



J2 W{H', H)P{H\ t)-J2 W{H, H')P{H, t), 



tion of the height profile h(x,t) is a noisy version of the 
Kuramoto-Sivashinsky (KS) equation and reads 



dh 

9^ = "° 



u\V'^h- K{V^fh+^{VhY +ri{x,t), (2) 



H' 



H' 



(1) 



where vq, v, k, and A are constants, and ri{x,t) is a 
Gaussian white noise. The KS equation (obtained when 
1] — 0) is often considered a paradigm of spatiotem- 
poral chaos and appears frequently in studies of pattern 
formation [ pd] |. An initially fiat one-dimensional inter- 
face described by the KS equation evolves in time from 
an almost periodic morphology at early times to a rough 
surface at late times [p^-p^ described by the Kardar- 
Parisi-Zhang (KPZ) equation 

The model for ion sputtering in Ref. |^] shows an ini- 
tial periodic morphology and a late-time KPZ scaling 
regime which are similar to the behavior of the KS equa- 
tion. However, the KS equation is a deterministic equa- 
tion whereas the erosion model is inherently stochastic. 
In |9| the early and late time dynamics of the erosion 
model were numerically studied with the conclusion that 
they are the same as those obtained from the noisy KS 
equation p^ . Here, we confirm this result by showing 
analytically that the noisy KS equation yields the con- 
tinuum description of the erosion model. By using the 
master equation approach, we determine the contribu- 
tions of the intrinsic noise to the evolution equation for 
this simple pattern forming model. 



where W{H,H') denotes the transition rate per unit 
time from configuration H to H' , and the sums are over 
all configurations H'. One can study the associated 
Fokker-Planck equation and therefrom derive a contin- 
uum Langevin equation for the surface height provided 
the fluctuations in the system are not too large, and pro- 
vided the continuum limit can be justified An im- 
portant feature of the master equation approach is that 
it provides a relation between the fiuctuations (noise) in 
the system and the parameters in the model. 

In the present paper, we follow the master equation 
approach in order to derive the continuum equation for 
a simple pattern forming model recently introduced to 
study the morphology of surfaces eroded by ion sput- 
tering 0. The equation which we obtain for the evolu- 



II. EROSION MODEL 

The model introduced in Ref. 1^ considers an interface 
in H- 1 dimensions described by the height variable hi (t) , 
i = 1, . . . ,L, where L is the system size. The sites be- 
low the interface are occupied with particles, whereas the 
sites above are empty. The lattice unit distance along the 
horizontal direction is denoted by a. The derivation of 
the noisy KS equation can be generalized to d+1 dimen- 
sions but in the following we discuss the 1-1-1 dimensional 
case. 

The two basic physical mechanisms acting on the sur- 
face are erosion and diffusion. In the model, a particle 
at the interface is chosen randomly. Then, the particle is 
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subjected to an erosion rule with probability /, and to a 
diffusion rule with probability 1 — /. 

The erosion rule is as follows: the particle at hi is 
eroded (the corresponding lattice position emptied) with 
probability PeYi The quantity is computed as 

1/7 times the number of occupied sites in a box of size 
3x3 centered at the chosen site, i.e., Pe is the num- 
ber of occupied nearest and next-nearest neighbors of the 
particle at hi, and takes the values 1/7, 2/7, . . . , 1. The 
definition of Pe accounts for the unstable erosion mech- 
anism which exists in the physical systems due to the 
finite penetration depth of the bombarding ions into the 
eroded substrate |Q. The efficiency of the sputtering 
process is measured by the sputtering yield [n9| 



Yi = Y{ip,) ^yo + yi ifl + y2 iff, 



(3) 



with the local slope (pi — ta.r^^^{{hi-^-l — /ii_i)/2a). In 
1^ yo, yi, and 7/2 were chosen such that Yi{0) — 0.5, 
Yi(TT/2) = 0, and with a maximum value Yi{ipm) = 1 for 
fm — 1 rad. In general one would merely require that 
yo >0,yi> 0, and < |l9). 

The surface diffusion rule moves a particle at the inter- 
face from the top of column k to the top of a randomly 
chosen nearest-neighbor column with the rates 



1 



l + exp(/3AHfe^fe±i)' 



(4) 



where the energy is H = {J/b^)Etli'^^ " ^'+1)' @' 
and ATi.k^k±i is the energy difference between the final 
and initial configurations; J is a coupling constant, b is 
the unit lattice spacing in the height direction, and (3 is 
the inverse temperature. 



the ratio if("+i)/ii:(") will typically be of order 0{l/n) 
|p2| . If the intrinsic fluctuations are sufficiently small 
then only the if term in Eq. (^) will survive in the 
limit — > 00 ||2^ . However, in our case the intrinsic 
fluctuations are relevant. Then, in the limit Q 00 the 
Kramers-Moyal expansion reduces to the Fokker-Planck 
equation 



dP 
'dt 



d 

dhi 



(7) 



which is obtained from Eq. (||) by keeping only the K^^'^ 
and terms 

The Fokker-Planck equation can be used in the further 
analysis of the system. However, only average quantities 
such as, e.g., {hi{t)) = '}2,^hiP{H,t), can be calculated. 
Therefore it is generally more convenient to recast Eq. (|^) 
in an equivalent Langevin form for the height hi(t). In 
the Stratonovich interpretation, the Langevin equation 
associated to the Fokker-Planck equation is 



dh 
'dt 



Vi 



(8) 



where terms of order 0{l/il) have been neglected, cf. 
p3[ . In Eq. (^), hi is a continuous variable describing 
the dynamics of the fluctuations of the height configura- 
tions {hi{t), . . . , hL{t)}, as obtained in the macroscopic 
limit Q 00, cf. (see also [p|,pO[). The term rji is 
a Gaussian white noise with average value equal to zero, 
and variance 



(9) 



III. STOCHASTIC FORMALISM 

The interface configuration at a given time t is speci- 
fied by the set of column variables H — {hi, h2, ■ ■ ■ , h^}. 
The moments of the transition rates W{H, H') in Eq. (|l|) 
are defined as follows: 



(5a) 



H' 



K. 



= Ee^^' - - ho)W{H', H), (5b) 

H' 



H' \k=l 



The first step in order to derive the Langevin equation 
is to transform the master equation (|^) into a Kramers- 
Moyal partial differential equation (see, e.g., |2l|) 



dP^^ 
dt ~ 



i-iy 



n\ dhi^ . . . dhi 



{k^I,^p) , (6) 



where a sum over repeated indices is assumed. 

Next, one identifies a parameter such that Q —>■ 00 
denotes the macroscopic limit. For \h^ — hi\ ~ 0{l/fl), 



IV. DISCRETE LANGEVIN EQUATION 

We now apply the formalism described in the previous 
section to the erosion model in The transition rate 
for the erosion rule reads 

We{H, i/') = - E PeYkS{h'^, hk-b)Y[ S{hr,hj), (10) 

where r is the time scale. The erosion probability Pe can 
be expressed in different forms. We choose 



ao 



(11) 



with the discrete Laplacian V^/i^ = ft-^+i — 2hi + hi^i, 
and constants ao and ai. In Eq. ( pi] ) the Laplacian ac- 
counts for the physical mechanism of the box rule, and 
the 8 term accounts for the finite size of the box. Defin- 
ing 9{x) = 1 for a; > 0, and 6{x) = for x < 0, the Q 
term can be written 

e = - {h,^i -h,- 2b) -h,-b] 

- {h,+i - h, - 2b) [h,+i - h, - b] 

+ 9 (h, - h,^i ~ 3b) [h,-h,^i-2b] 

+ 6{h,- h,+i - 36) [h, - h,+i - 2b] . (12) 
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The 9 term in Eqs. (|T^) and (|T2|) becomes effectively zero 
in the hydrodynamic hmit of the model where the surface 
roughens in a way consistent with the KPZ universality 
class and the slopes are small along the interface. Also, 
from the numerical simulation it is known that the 
qualitative behavior of the model does not change if the 
size of the box is enlarged [E5|. The sputtering yield is 



obtained from Eq. (||) as 



2/2 - 2yi/3 



with the discrete gradient V/i^ — (/li+i — ft.i_i)/2. 
The transition rate for the diffusion rule reads 



(13) 



Wd{H, H') 



1-/ 
2t 



Y,HSihlhk~b)6iK+,,hk+i+b)+w^^,5ihlh, + b)S{h',+,,hk+i-b)] Yl 5ihr,h,), (14) 



j^k,k+l 



where are defined in Eq. (^. By expanding the transition rates w^, we find 

1 2l3Jq/b + E ^ (^''^'^ " ^"h^^^T ^ 



1 + 9 (1 + 9) 



n=2 



(15) 



where q = exp(6/9J), and C„(g) are numerical constants 
whose exact value will not be needed in what follows. 

Using that W{H,H') = We{H,H') + Wd{H,H'), and 
expression (^, we obtain to lowest order the following 
values for the transition moments: 



K, 



(W)Mv2(V^/.,0, (16a) 



p y — 



(2)_ {l-fy 



K ' = 



T{l + q) 



fb' 



(16b) 



Here, 



of the form (V^/i- 



26,, 



i,j + <5i-ij- Higher order terms 
V^/ii±i)^"+^ have been omitted (in 
section ^ we will argue that such terms are irrelevant 
in the renormalization group (RG) sense for the scaling 
properties of the interface). The quantity PeYi reads 



PeY. = ^ + "^^'h + |^(V/».)', (17) 



5yo 
7 



with additional terms of the form {Vhi)'^™'V'^hi. 

The different terms appearing in Eq. ( [l^ ) give rise to 
the following contributions in the Langevin equation (|^): 
an additive constant which corresponds to the average ve- 
locity of erosion for a flat interface, a negative Laplacian 
which reflects the unstable nature of the erosion rule, and 
a iyhiY' nonlinearity of the KPZ type describing lateral 
motion of the interface JTst ; the nonlinearity originates 
from the coefficient yi in expression (^) for the sputter- 
ing yield. In addition, the diffusion rule contributes to 
Eq. (H) with the linear term V^/ii (plus higher order non- 
linearities). As demonstrated in Ref. by numerical 
simulations and shown using a coarse-graining procedure, 
the V^/ii term originates from the conserved nature of the 
surface diffusion rule (||) [||j2^ . 



V. CONTINUOUS LANGEVIN EQUATION 

To obtain a continuous Langevin equation, we assume 
that the discrete function hi(t) can be replaced by a 



smooth function h{x,t) such that hi[t) — h{x — ia,t), 
and further that 



(18) 



A more rigorous way to arrive at the function h[x,t) 
would be through some coarse graining procedure which 
preserves the symmetry of the problem. However, as dis- 
cussed in (|,|2^ , this is a highly nontrivial task. Here we 
are concerned with the form of the relevant terms in the 
equation of motion for h(x, t), which are expected to have 
the same form as those obtained by the use of (|l8|). We 
will determine them as the leading terms in an expan- 
sion in which we take the parameters a and b to be small 
but nonzero 1^. Combining (|l|) with Eqs. (|), (|) we 
obtain the continuous Langevin equation for h{x,t): 



dh 

at 



= Vq 



K{V^fh+^{Vhf +'q{x,t) (19a) 



+ci(V/i)2v2/i + C2{'^^hfV^h + . . . , (19b) 

where now V = d/dx. The noise •q{x^t) has correlations 

(?7(x, t)r]{x',t')) = [2D - 2DdV^] d{x - x')d{t ~ t') 

+ [DiV^h + D2{Vh)'^] 6{x - x')6(t - t')- (20) 

The coefHcients in Eqs. (|l|)-(|20l) read 



"0 



5fyob 
7t ' 

lOfyib 

7t ' 
fyoaapb'^ 
7t ■■ 



D = 



Do 



fyoapb 
7t ' 

Ut ' 

5fyiab^ 
7t ' 



(1 - f)pjqa^ 



D, = 



T{l + q? ' 
(1 - f)a%^ 



2T(l + g)2 ' 
(21) 



The terms in Eq. ( |l9a|) constitute the noisy Kuramoto- 
Sivashinsky equation (^ , where one notes the presence of 
a negative Laplacian and a KPZ nonlinearity. The pres- 
ence of such a nonlinearity is known to determine the 
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scaling behavior of the noisy KS equation as shown by 
an RG analysis ||2^. Using the values of the roughness 
and dynamic exponents at th e KPZ fixed point ITq], the 
terms appearing in Eq. ( 19b ) can be shown to be irrel- 
evant. A similar RG argument shows the irrelevance of 
the conserved term Dd in the noise correlator ( po| ) in the 
presence of the shot noise D, as well as the irrelevance of 
the multiplicative contributions Di and D2- 



properties as the Kardar-Parisi-Zhang equation and this 
result is in nice agreement with the simulations of the 
erosion model performed in Ref. The master equa- 
tion approach has allowed us to derive the form of the 
stochastic noise in the Langevin equation for a simple 
pattern forming system. 
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In summary, we have derived the continuum equation 
for a microscopic model for ion sputtering. The resulting 
equation is a noisy version of the Kuramoto-Sivashinsky 
equation. Such an equation exhibits the same scaling 
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